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Abstract 

We develop a model of magnetoresistance oscillations induced by the Hall 
field in order to study the temperature dependence observed in recent exper- 
iments. The model is based on the solution of the von Neumann equation 
incorporating the exact dynamics of two-dimensional damped electrons in 
the presence of arbitrarily strong magnetic and dc electric fields, while the 
effects of randomly distributed neutral and charged impurities are perturba- 
tively added. Both the effects of elastic impurity scattering as well as those 
related to inelastic processes play an important role. The theoretical predic- 
tions correctly reproduce the main experimental features provided that the 
inelastic scattering rate obeys a T 2 temperature dependence, consistent with 
electron-electron interaction effects. 



1. Introduction 

Magnetoresistance oscillations of two-dimensional electron systems (2DES), 
besides that Shubnikov-de Haas oscillations, have been a considerably active 
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Figure 1: (color online). Single quantum scattering rate l/r g vs T 2 



topic of research over the past decade. Microwave-induced resistance oscil- 
lations (MIRO) were registered 0, 0, 0] in high mobility semiconductor 
2DES samples subjected to microwave radiation and low magnetic fields. 
Similarly Hall field-induced resistance oscillations (HIRO) are observed in 
high mobility samples in the presence of an strong longitudinal dc-current 
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Even though the mechanisms that produce MIRO and 
HIRO are different, both display periodic oscillations in 1/B that originates 
in the commensurability of the cyclotron frequency u c with a characteristic 
parameter of the system. Remarkably, in sufficiently clean samples the MIRO 
AND HIRO evolve into zero resistance or zero differential resistance states, 
respectively 0, [loj]. Some models have been proposed to explain the main 
experimental features of HIRO. On one hand the displacement model relays 
on the impurity induced electron transitions between Landau levels (mod- 
ified b y th e effects of the microwave radiation or the strong dc-excitation) 
@, EH 12 1 on the other, the inelastic model depends on the formation of a 
non-equilibrium distribution function induced by the Hall field flif . 

Recently there has been considerable interest in the origin of the temper- 
ature dependence of MIRO 14|, |l5|, [l6| and HIRO 17|. In both cases the 



oscillations are best observed at T ~ 1 K and are smeared when T is of the 
order of a few Kelvin. These experiments show that the amplitude of the 
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Figure 2: (color online). Longitudinal differential resistivity function of the 

inverse of the magnetic field 1/B for temperatures from T = IK (thick [blue] line) to 
T = 6K (thick [blue] line) in steps of 0.5K (thin [gray and red] lines). 
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Figure 3: (color online). Normalized amplitude of the HIRO oscillations Ar/po as a 
function of the inverse magnetic field 1/B for the temperatures T = 1, 3, 4, 4.5, 5, 5.5, 6K 
for cti = 25. The points are obtained from the oscillation maxima of the plots in Fig. 2,, 
whereas the lines are the predicted linear behavior of the Dingle factor, Eq. (|22l . 
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Figure 4: (color online). The normalized HIRO amplitude (normalized to the offset of the 
Po) Ar/po as a function of 1/lo c t s . The circles indicate ai = 15 and a,; = 25, temperatures 
ranging from IK to 6K in 0.5K intervals and magnetic fields between 0.02 and 0.2T. In 
this plot all data points obtained from our model for different temperatures, magnetic field 
and a values fall on the same universal line (solid line) prescribed by Eq. ([21) . 

oscillations scale as exp(—aT 2 /B). In this paper we present a model that 
correctly reproduces the temperature dependence of HIRO. Both the effects 
of elastic impurity scattering as well as those related to inelastic processes 
due to electron-electron scattering play an important role. The model is 
based on the exact solution of the von Neumann equation for a 2D damped 
electron gas subjected to arbitrarily strong magnetic and dc electric fields. 
In addition the model incorporates the weak effects of randomly distributed 
neutral and charged impurities through second order time dependent pertur- 
bation theory. This procedures yields a Kubo-like formula for the electric 
density current that incorporates the non-linear dependence on the electric 
field. Both inter- and intra-Landau level transitions contribute to the density 
current. 



2. Model 



Our starting point is the Hamiltonian for an electron in the effective 
mass approximation in two dimensions, subject to a uniform perpendicular 
magnetic field B = Bk, an in-plane electric field E = E x i + E y j, and 
the impurity scattering potential V(r). Therefore the electron dynamics is 
governed by the total Hamiltonian H = H e + V, with 

H e = Ho + eE-x, (1) 

here H = II 2 /2m, m is the effective electron mass, e is the electron's charge, 
II = p + eA is the velocity operator and the vector potential for an homoge- 
neous magnetic field in the symmetric gauge is given as A = (—By,Bx) /2. 
We consider the impurity scattering potential to be composed of a com- 
bination of short-range uncorrelated neutral delta scatterers and charged 
scatterers 

V (r) = e^l £ / T^Vn (q) exp [iq ■ (r - r*)] 

+ e_r|t| E J J^f V < (?) ex P [«'« • ( r " p i)l ' ( 2 ) 

where r« is the position of the zth neutral scatterer, is the position of the 
jth charged impurity and N n and N c are the number of neutral and charged 
impurities. The explicit form of the Fourier coefficients V n (q) and V c (q) 
depends on the nature of the scatterers. Neutral scatterers are characterized 



by a constant (la, ll9|, ll2 



Tr , , 2irh 2 a 

V n (q) = , 3 

m 

where a is a dimensionless parameter. The explicit calculations only allows 
us to calculate the product a 2 N n , hence, without lost of generality, it is set as 
a = 1. The potential for charged impurities is modified by screening effects 
and takes the form 

fe2 -qd 

Qtf 

where d is the thickness of the doped layer, e& is the relative permittivity of 
the surrounding media, and the Thomas-Fermi wave number is given by 

e 2 m 
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The adiabatic switching of the impurity potential at the initial time t — > — oo 
is controlled by the factor exp (— T\t\). 

Now we turn to the calculation of the density current. It is computed as 
the time and thermal average of the velocity operator 



TS 



T/2 

dm [ P (t) n] 

■T/2 



(6) 



T->oo 



where S is the sample surface, and p (t) is the density matrix operator ob 
tained from the solution of the von Neumman's equation ihdp/dt = \H, p] 



This equation is solved by means of a series of unitary transformations [23 
The magnetic and electric field part of the Hamiltonian are treated exactly 
by means of the first unitary transformations, that is expressed in terms of 
functions that solve the corresponding classical equations. Whereas the im- 
purity potential is incorporated through second order perturbation theory. It 
is assumed that, in the absence of the impurity, potential the density matrix 
reduces to the equilibrium density matrix given by po = f(H e ), with f(E) 
the Fermi distribution function. 

The computation of Eq. fl2]) yields an explicit expression for the current 
density. It splits into the Drude and an impurity induced contribution: 

J = J D + J im P . (7) 

The Drude contribution is given by J D = n e e Vd, where n e is the electron 
density and the drift velocity is given by 

_er™ E-u c r™kxE 



c'tr 



in 
tr 



The transport inelastic scattering time r t l " arises from a damping term v/r t 
that is added to the classical equation of motion in order to incorporate 
dissipative effects. The components of impurity induced density current can 
be expressed as 

Jim = e -wllf d ^ [W - h,-,/*] K<w + ( 9 ) 

mm' 

where n n = N n /S and n c = N c j S are the density of neutral and charged 
impurities; / is the Fermi distribution function evaluated at the tilted LL 
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energies £^ q = huj c (n + 1/2) + hu q where u q = q ■ v d . The function G w ' is 
given by 

<W = \Vn,c (g)| lA^ (%)] A2 - ~ 4cj2r2 , (10) 

where q = (q x ,Qy) is the transfer momentum, and q = (q y ,—q x ) its dual, 
whereas 

A^ = [co q + co c ^-fi')} 2 -tu 2 c +T 2 . (11) 
Finally the matrix elements D^ v are given by 




z 



L u> M (\ Z q\ ) ^ > 



IV CO = exp -i^J_ y " v /« " - ' ' (12 



M-' l {\z q \ 2 ),IM<ljt, 



where z g = /_b(<7x — iqy)/V2, and L^ -1 ' denotes the associated Laguerre poly- 
nomial. 

Though J D depends linearly on the electric field, J imp introduces a non- 
linear contribution through u q . From the Lorentzian form of the G n £, func- 
it can be seen that the T parameter effectively controls the 



tion in Eq. ([TO 

LL broadening |24j|. We thus introduce disorder broadening effects through 



T, moreover, in order to take into account the LL width dependence on the 



magnetic field henceforth we shall consider [25l. |26 



r = J— (is) 



7TT q 



as expected we assume that the LL width depends on the quantum scattering 
times r q . 

The differential conductivity tensor is calculated from = dJi/dEj and 
the differential resistivity tensor is obtained from the inverse of the conduc- 
tivity: that is Ty = cr^j . 

HIRO experiments are carried out at a fixed longitudinal current density 
Jdc and zero transverse current. The longitudinal electric field E x and the 
Hall field E y are thus computed from the following implicit equations 

J x (E x , E y ) = J dc , J y (E x , Ey) = 0, (14) 

where J x and J y are given by Eqs. flTl fTT]) . Given the complexity and the 
nonlinear character of the current density expressions, the solution of Eq. 
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(TTTj) is computed by a recursive numerical iteration, until the accuracy goal 
J y /J x «lx 10~ 15 is reached [12 . 

In the present model the current density Eq. (J7J) depends directly on the 
transport inelastic scattering time Eq. (jSJ), and the quantum scattering time 
Eq. (fl3|) . On the other hand, in order to make contact with experimental 
results we need to specify the mobility, that is given in terms of the total 
transport time \x = er tr /m. Using the Matthesiessen's rule it is possible 
to incorporate both inelastic and impurity contributions to the transport 
scattering time as 

1 1 

(15) 



1 

Ttr 



T, 



tr 



T, 



imp 



tr 



The transport as well as the quantum scattering times associated to the 



impurity potentials are given by[18l. 120 



_imp 



m 



2kF , \V n (q)\ 
■ / dq 





1 _ f-2_)2 
1 V 2k F > 



: 2kF , \V c (q)\ 
+ n c dq 

lo Jl 



Q \2 
2k F > 



(16) 



and 



m 



^ 2, \Vn{q)\ 

n n I q dq 

o 



1 _ f_2_)2 

1 V 2k F > 



lo Jl 



_Q_\2 

2k F > 



(17) 



where n n and n c are the neutral and charged impurity densities respectively, 
and kF is the Fermi momentum. For the quantum inelastic time we assume 
that it arises from electron-electron scattering, that is known to be well 



estimated by the following expression [2ll . |22 

2 



1 



Ef ( k B T 
h { E F 



+ 1 



(18) 
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For typical experimental conditions, the following condition: In (if~r) 
In ^2^f J <C 1 holds; thus l/r* n has mainly a T 2 behaviour 



+ 



1 E F fk B T x 2 



(19) 



We have explicit expressions for the quantum impurity Eq. ( TT6|) and quantum 
inelastic Eq. ( TT9|) scattering times; hence it is straightforward to determine 
the total quantum scattering rate using 

iii 

(20) 



_ imp --in 

The determination of the transport inelastic time is more involved given 



the complex structure of the collision integral [22|]. Thus, for simplicity, we 
shall assume a simple relation relation between the transport and quantum 
inelastic scattering times: t™ = oiiT™ where a{ is a constant. Now we can 
use Eq. (fT5j) to determine the total transport time. 



3. Results and discussion 

We consider a sample of electron surface density and mobility of n e = 
3.7 x 10 11 cm~ 2 and /i = 1 x 10 7 cm 2 /Vs respectively, a doped layer of d = 
20nm with n c /n n = 300 and n n = 1.7 x 10 7 cm -2 that leads to a ratio between 
the quantum and transport impurity scattering times of rl™ p /T l q mp ~ 11 
according to Eqs. f[T6]) and fTT7|) . The longitudinal current / is set to 80/iA 
and for the sample width we assume w = 100/im; hence Jd c = 0.8 A/m. 

We have a formalism in which the nonlinear resistivity effects, as well 
as the relaxation rates can be consistently calculated, once the neutral and 
charged scattering potentials have been specified. The quantum scattering 
rate l/r q as a function of T 2 is presented in Fig. [TJ As the temperature 
increases, we observe the linear behavior related to the inelastic scattering 
contribution Eq. (I19|) . On the other hand, the low temperature limit is de- 
termined by the impurity scattering contribution in Eq.f ll6p . 

In Fig. [2] we present the longitudinal differential resistivity r xx as a func- 
tion of the inverse magnetic field for temperature values from IK to 6K in 
0.5K increments. We observe clear differential magnetoresistance oscillations 
periodic in 1/B, we see up to the 6th cyclotron resonance. The HIRO origi- 
nate from the inter-Landau level transitions induced by impurity scattering 
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and are governed by the ratio of the Hall to the cyclotron frequencies. We 
observe that the oscillation peaks appear near integer values of the dimension- 
less parameter e = ojh/u c . Here Hujh = gEh(2R c ) is the energy associated 
with the Hall voltage drop across the cyclotron diameter R c = vp/u c and 
Eh = JdcB/en e is the Hall field. The oscillatory periodic pattern is clearly 
temperature independent, however we observe that the oscillation ampli- 
tude gradually decay with increasing temperature, and almost disappear at 
T PS 6K. To further characterize the HIRO temperature dependance as well 
as its physical origin, we present in Fig. [3] the normalized HIRO amplitude 
Ar/po = ( r xx — ro) / po corresponding to the same data shown in Fig. flTJ 
and arrange it in a semilog plot as a function of the inverse of the magnetic 
field for temperatures ranging from 2K up to 6K in IK steps in. Here po is 
Drude resistivity and r$ is the offset of the HIRO. It is interesting to observe 
that the oscillations amplitude clearly display an exponential decay. One 
also observe that the exponent grows with increasing temperature as T 2 and 
decreases with the magnetic filed as 1/5. The results of our detailed calcula- 
tion reproduce the most important features of recent experiments and is also 
in g ood agreement with the approximated expression proposed in reference 
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r X x = r + p cos 27r — (21) 

where T n is the time describing electron backscattering off impurities and S 
is the Dingle factor 

6 = exp(-—), (22) 

given in terms of the total quantum scattering time Eq. (l20p . 

Our complete calculation depends on the total quantum scattering rate 
through the LL width Eq. (TT3]) and also on the transport inelastic rate 
that appears on the term Drude 3 D = n e evd with the drift velocity given 
by Eq. (JSJ) . Consequently the differential resistivity is expected to depend 
on the temperature, magnetic field and the parameter ccj that relates the 
transport and quantum inelastic scattering times: t™ = otiT l q n where CKj. 
However Eqs. (12 1 |I22[) suggest a simple relation for differential resistance on 
terms of the Dingle factor. In Fig. H] we present the results of our model 
for the normalized HIRO amplitude Ar / p normalized to the offset of the po 
HIRO oscillations as a function of 1/ uj c t q where the different symbols (colors) 
indicate «j = 15 and on = 25, temperatures ranging from IK to 6K in 0.5K 
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intervals and magnetic fields between 0.02 and 0.2T. In this plot all data 
points obtained for different temperatures, magnetic field and a, values fall 
on the same universal line prescribed by Eqs. (j21|22p . Hence it is concluded 
that the HIRO temperature dependance originates from electron-electron 
contribution to the quantum scattering rate. 

4. Conclusions 

We have presented a model to explain the HIRO properties of high mobil- 
ity 2DES, including its temperature dependence. The effects of the electric 
and magnetic part of the Hamiltonian were treated exactly, whereas the im- 
purity potential is included by means of a second order perturbation calcula- 
tion. Our model contains the effects of neutral delta scatterers and charged 
impurities. We introduced inelastic effects through the damping that ap- 
pears in the solution of the classical equation of motion that leads to the 
Drude contribution Eq.flH]). The T 2 dependence of r* n is required in order to 
reproduce the correct temperature dependence of the HIRO amplitude, and 
is consistent with electron-electron inelastic processes. 

We acknowledge support from UAM-A CBI project 2232204. A. Kunold 
wishes to thank INSA- Toulouse for the use of their cluster. 
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